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Abstract 

The generating function for the (colored) Motzkin numbers of higher rank introduced recently is 
discussed. Considering the special case of rank one yields the corresponding results for the conventional 
colored Motzkin numbers for which in addition a recursion relation is given. Some explicit expressions 
are given for the higher rank case in the first few instances. 



1 Introduction 

The classical Motzkin numbers (A001006 in P^) count the numbers of Motzkin paths (and are also related 
to many other combinatorial objects, see Stanley [2]). Let us recall the definition of Motzkin paths. We 
consider in the Cartesian plane Z x Z those lattice paths starting from (0,0) that use the steps {U,L,D}, 
where U — (1,1) is an up-step, L — (1,0) a level-step and D = (1,-1) a down-step. Let M{n,k) denote 
the set of paths beginning in (0,0) and ending in (n,k) that never go below the x-axis. Paths in M{n,0) 
are called Motzkin paths and m„ :— |M(n,0)| is called n-th Motzkin number. Sulanke showed ^ that the 
Motzkin numbers satisfy the recursion relation 

{n + 2)m„ = (2n + l)m„_i + 3{n - l)m„_2 (1) 

and it is a classical fact (see Stanley j2j) that their generating function is given by 



l^rrinx = . (2) 

n>0 

Those Motzkin paths which have no level-steps are called Dyck paths and are enumerated by Catalan numbers 
(A000108 in [T]), see Stanley 2]. In recent times the above situation has been generalized by introducing 
colorings of the paths. For example, the k-colored Motzkin paths have horizontal steps colored by k colors (see 
[HIS] and the references given therein). More generally, Woan introduced '6'i'7^ colors for each type of step. 
Let us denote by u the number of colors for an up-step U , by I the number of colors for a level-step L and 
by d the number of colors for a down-step D. (Note that if we normalize the weights as u + l + d = 1 we can 
view the paths as discrete random walks.) One can then introduce the set M*^"''''^)(7i, 0) of {u, I, d)- colored 



^AU correspondence should be directed to this author. 
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Motzkin paths and the corresponding {u,l,d)-Motzkin numbers m^"'''''^ :— |M("'''''^(n, 0)|. Woan has given 
[6] a combinatorial proof that the (1, ^, c?)-Motzkin numbers satisfy the recursion relation 

{n + 2)m(i^'-^) = l{2n + l)mi':!f + {M P){n l^^^^f. (3) 

Choosing I = 1 and d ~ 1 yields the recursion relation ([T]) of the conventional Motzkin numbers to„ = 
^(1,1, i)_ Note that choosing {u,l,d) = (1, fc, 1) corresponds to the fc-colored Motzkin paths. Defining 
mk^n ■= |M(^''''^)(n, 0)1, one obtains from ([3]) the recursion relation (n + 2)mk,n = k{2n + l)mfc^„_i + (4 — 
k'^){n — l)mfc_„_2 for the number of fc-colored Motzkin paths. Sapounakis and Tsikouras derived Jt^ the 
following generating function for mfc.„: 



„ 1 - fc.x - - kxY - 4a;2 
E = ■ (4) 

ri>0 

For fc = 1 this identity reduces to ^ for the conventional Motzkin numbers m„ = mi_„. One of the 
present authors suggested [S] (as "Problem 1") that it would be interesting to find the recursion relation 
and generating function for the general (u, /, (i)-Motzkin numbers ml"''''^^ We will prove in Theorem 12.11 
thet ml"'''''' = mn''''"'^\ yielding the desired recursion relation. Furthermore, a generating function and an 
explicit expression is derived for ml"''''''. In |8] it was furthermore suggested to find the recursion relation 
and the generating function for the (colored) Motzkin numbers of higher rank ( "Problem 2" ) . These numbers 
were introduced by Schork ^ in the context of "duality triads of higher rank" (and have also been considered 
before as "excursions", see, e.g., [IH])- In view of this connection Schork conjectured [Hj that the Motzkin 
numbers of rank r satisfy a recursion relation of order 2r-|- 1. Very recently, Prodinger was the first to observe 
pT] that this conjecture does not hold by discussing explicitly the case r = 2. In particular, already for the 
first nontrivial cases r = 2, 3 the relations involved become very cumbersome. We will describe in Theorem 
13.21 the generating function for the Motzkin numbers of higher rank and discuss then several particular cases 
explicitly. 

2 Recursion relation and generating function for the general (col- 
ored) Motzkin numbers 

Theorem 2.1. The general {u,l,d)-Motzkin numbers satisfy the recursion relation 

{n + 2)mt'''''> = l{2n + l)m^^l'f + {And - f)in - l)m^^lY'^ . (5) 
A generating function is given by 

V _ l-lx- ^jl-lx)^ ^4udx^ 

^ " 2udx^ ' ^ ' 

n>0 

implying the explicit expression 
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Proof. Let us prove that mn"'''''^ = run''''"'''^ for all n > 0. In order to see that, let U be any up-step in a 
Motzkin path; we call {U, D) a pair if the down-step D is the first down-step on the right-hand side of U 
which has the same height as U. The set of pairs of a Motzkin path is uniquely determined and for each 
pair exist ud possible combinations of colorings. The same number of possible colorings result if the up-steps 
u are always colored white (i.e., u = 1) and if each down-step can be colored by ud colors (which we call 
"alternative colors" to distinguish them from the original colors). Thus, given a (w, /, (i)-Motzkin path we 
may replace the colors {uj, dj) (with I < Uj < u and 1 < d-, < d) for the j-th pair ([/, D) by the combination 
of alternative colors (1, Cj) where cj is the {[uj — l]d + dj)-th. alternative color (with 1 < Cj < ud). Replacing 
the colors of all pairs in this fashion by the alternative colors yields a (1, 1, M(i)-Motzkin path. Thus, we have 
constructed a bijection between the set of (u, (i)-Motzkin paths of length n and set of (1, ^, ud)-Motzkin 
paths of length n, thereby showing that rri'n'^''^'^ — rrin'^'^'^K From this identity and ([3]) we immediately 
obtain ([5]). An equation for the generating function M(-„ ; ,^-)(a;) := X]ri>o obtained from the 
"first return decomposition" of a nonempty (u, (i)-Motzkin path M: either M — LM' or M — UM' DM" , 
where M',M" are (u, (i)-Motzkin paths. The two possibilities give the contributions lxM(^u i ,i){x) and 
udx^(M(„^;^(i)(x))^. Hence, M(„_;^d)(a^) satisfies 



yielding 



l~lx- J{l-lxy-4udx'^ 



which is the asserted equation ^ for the generating function. Note that this may also be written as 



where 

is the generating function for the Catalan numbers (see Stanley [5]). Thus, 

1 f2j\ uH'x"^^ 

Recalling M(„ ; (^^(a:) = J2n 'tiL"''''^^^", a comparison of coefficients shows that the number of (u, I, d)-Motzkin 
paths of length n is given by ([7]) . □ 

3 Generating function for the general (colored) Motzkin numbers 
of higher rank 

We will now generalize the situation considered in the previous section to the case of higher rank. One of the 
present authors discussed [S] in the context of duality triads of higher rank (where one considers recursion 
relations of higher order, or equivalently, orthogonal matrix polynomials '12|) why it is interesting to consider 
the situation where the steps of the paths can go up or down more than one unit. The maximum number of 
units which a single step can go up or down will be called the rank. More precisely, let r > 1 be a natural 
number. The set of admissable steps consists of: 
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1. r types of up-steps Uj — (1, j) with weights uj for 1 < j < r. 

2. A level-step L = (1, 0) with weight I. 

3. r types of down-steps Dj = (1, —j) with weights dj for 1 < j < r. 

In the following we write (u, I, d) :— {ur, . . . , mi, /, di, . . . , d^) for the vector of weights. 

Definition 3.1. [9 The set M^^'^''^\n, 0) of (u, I, d)-colored Motzkin paths of rank r and length n is the set 
of paths which start in (0, 0), end in (n, 0), have only admissable steps and are never below the a;-axis. The 
corresponding number of paths, ml"''"'*'' :— |M'^"''^'^^(n, 0)|, will be called (u, Z, d)-Motzkin number of rank r. 

Motzking paths of higher rank were considered in the literature already before [9] under the name 
"excursions", see, e.g., (TU]. In [5] it is discussed how one may associate with each Motzkin path of rank r 
and length n a conventional Motzkin path of length rn in a straightforward fashion. However, it was also 
discussed that this association is no bijection and that the case of higher rank is more subtle. In the following 
we derive an equation that the generating function for the Motzkin numbers of higher rank satisfies. In order 
to do that we need the following notations. We denote the set of all (u, I, d)-Motzkin paths of rank r that 
start at height s and end at height t (and have only admissable steps and are never below the x-axis) by 
As,t, and we denote the subset of paths of length n by As,t{n)- Define 

As,t = As^t{x) := \As,t{n)\x". 

n>0 

We extend this notation by defining Ag^t = for all s < or t < 0. For s, i > we denote the subset of 
paths in As.t that never touch the a;-axis by Al j. In the following we make several times use of the relation 

■^*s,t - As-i,t-i, or, \Al^t{n)\ = \As-i,t-i{n)\, implying 

Al,{x) I-^mWI^" = A-i,t-i{x)- (8) 

n>0 

Theorem 3.2. The generating function Aq^ = M(^j^ i^^j{x) — X]n>o f'^''" general (u, Z,d)- 
Motzkin numbers of rank r satisfies 

r r 

Ao,Q = 1 + IxAofl + x'^Aofl Y Y ^*pcig^P-i,g-i > (9) 

p=l g=l 

where the generating functions Aij with 1 < i < r — 1 and < j < r — 1 satisfy 

r 

AiJ — Ai—lJ — l + xAqj ^ ^ dqAi—l^q—l, 
q=l 

and for all I < j < r — 1, 

r 

Aoj = xAa^o ^ UpAp_ij_i. 
p=i 
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Figure 1: First return decomposition of a Motzkin path of rank r. 



Proof. From the definitions, a nonempty (u, Z, d)-Motzkin path M can start either by a level step L or 
an up-step Up with 1 < p < In the case the path starts by an up-step Up, we use the "first return 
decomposition'" of A/, i.e., we write M as M = UpM' DqM" where M" is an arbitrary (u, Z, d)-Motzkin 
path and M' e Ap^q such that M' does not touch the height zero, see Figure [TJ Thus, M' e Ap ^ and the 
generating function Aq.o satisfies 

r r 

^0,0 = 1 + IxAqa + a;^Ao,o ^ ^ UpdqAp^i^q^i 

p=l q=l 

where we have used (|S]). This shows Now, let us write an equation for the generating function Ai^j 
(with 1 < i < r — 1 and < j < r — 1) for the number of (u, Z, d)-Motzkin paths Qij = UiPijDj where 
Pij G Aij. In order to do that, we use the first return decomposition of such paths: either Qij — UiP^jDj 




Figure 2: The case 1 < i < r — 1 and < j < r — 1. 



such that Plj does not touch the height zero (i.e., P/^ G A*j), or Qij = UiP^qDqP^jDj such that Dq is the 
first down-step that touches the height zero (thus, Plq does not touch the height zero, i.e., P/^ G -^*g) 
described in Figure [2l Using ([8]), the generating function Ai^j thus satisfies 

r 

Aij = + xAo^j ^ dgAj_i_q_i, 

q=l 

as claimed. 

In order to write an equation for the generating function ^o.j (with < j < — 1)) let us consider the last 
up-step Up from height zero (there must exist at least one such up-step). In this case each (u, Z, d)-Motzkin 
path can be decomposed as P'UpP^jDj, where P' is an arbitrary (u, I, d)-Motzkin path and P^j E Ap,j such 




Figure 3: The case i — and 1 < j < i — 1. 
that it does not touch the height zero (i.e., Ppj £ Ap j), as described in Figure O Thus, using again ([5]), the 
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generating function Aq^ satisfies 

r 

p=i 

which completes the proof. □ 



Example 3.3. Let us consider as an example the case r = 1. It follows from Theorem \8.S\ that 

= 1 + (Lt - 1)^0.0 + uidix^Al^^, (10) 

yielding 



Ix — — Ix)^ — Auidia 



2uidix 
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as described in the proof of Theorem \2.1\ This implies - according to ^ - that 

7/2 



j=0 

Example 3.4. As another example, it follows from Theorem \3.Si that one has for r — 2 that 

Ao,o = 1 + IxAq^o + Ao.o{uidiAofi + uiii2^o,i + "2^1^1,0 + "2^2^1,1), 
^0,1 = a;^o,o(Mi^o,o + ■"2^1,0), 
Ai^o = xAo^o{diAofi + d2Ao^i), 

Ais = ^0,0 + a;Ao,i(diAo,o + ^2^0,1)- 
Solving the above system of equations we obtain that 

= 1 + {Ix — l)A()^o — x'^{d2U2 — diUijA^f^ + x'^{xu\d2 + xu2d\ — 2xlu2d2 + 2^2^2)^0,0 
-U2d2x'^{d2U2 - diWi)Ao Q + x'^u^dKlx - 1)AI q + x^u^d^AQ ^. 

Note that, if we set U2 = d2 = in the above expression then we get (jlOp . It is also interesting to consider the 
case (u, I, d) = (1, 1, 1) of non-colored Motzkin paths of rank 2, representing the most natural generalization 
of the conventional Motzkin paths. The above equation reduces in this case to 

= l + {x-l)Ao^o + 2x^Alo + x^{x-l)Al„+x^Alo 

= (1 + .tAo.o)'(1 -{x + l)Ao,o + xix + 2)Alo - x^{x + l)^^ q + x^^^o)- 

Setting Aq^ = M(^-i,i.i){x) and using that Aq^ is a formal power series (thus, Ao,o 7^ ~- ), this is equivalent 
to ' 

= 1 - (x + l)M(i,i,i)(x) + x{x + 2)M(2, 1 i)(x) - + 1)M(\ ^ ^^(x) + x^M^^ ^ i)(x), 

as described already in \1C\ To obtain a recursion relation for m„ = nrih''^''^\ one can use the MAPLE 
program package gfun written by Salvy et al. \13f^ . Prodinger has done this 111] to obtain 

625(n + 3)(n + 2)(n + l)m„ - 125(n + 3)(n + 2)(7n + 27)m„+i - 50(n + ?,){bn^ + 24n + 23)m„+2 
+(41890 + 30860n + 754071^ + 610n^)m„+3 - (6844 + 5151n + 121471^ + 91n^)m„+4 
-(n + 7)(23n^ + 301n + 976)m„+5 + 2(2n + 13)(n + 8)(n + 7)m„+6 = 
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n = l 


2 


3 


4 


5 


6 


7 


8 


9 


10 


(1.1,1) 

mil 


1 


3 


9 


32 


120 


473 


1925 


8034 


34188 


147787 



Table 1: The first few values of the non-colored Motzkin numbers of rank two. 



and mentions that Salvy has informed him that this recursion of order 6 is minimal. Thus, the Motzkin 
numbers of rank r = 2 satisfy a 7-term recursion relation and not a 2 • 2 + 1 = 5-term relation as the 
conjecture of Schork 181 implies. Thus, the conjecture does not hold! This was first observed by Prodinger 
\1 1]/ . The first few values of mn = m^'^''''^ are given in TaWe [TJ- this sequence is sequence AIO4I84 in '(Jj. 

The above example shows that the general case (where the weights are not restricted to 1) seems to be 
extremely complicated (see Theorem l3.2[) . Thus, from now on let us consider the case where all weights are 
equal to 1, i.e., u.j, = I — di — 1 for 1 < i < r. In this case we denote the set of paths by B^^t or Bs^t{n) and 
the generating function by Bi j (instead of A^j). Then one has 

Bs,t = Bt,s (11) 

which can be easily seen since to each path P G Bs.t{n) one can associate a path P' e Bt.s{n) by traversing 
P in opposite direction. Since this is clearly a bijection the above equation follows. Theorem 13.21 shows that 
one has to solve in the general case a system of equations for the unknowns Aij . Due to the symmetry 
([TT|) this reduces in the case where all weights are equal to 1 to a system of ' '"^2^'' equations in the ''^-''^"'"-^ 
unknowns Bij (where r — l>i>j>0). 

Example 3.5. Theorem \3.'2\ with pip gives for r — ?> the following set of ^ — 6 equations 





= 1 + xBq,o + x'^Bofl 


{Bo.o + 


2Bi,o + 2i?2,o + 2i?2,i + Si,i + -82^2), 


Bifi 


= xBofi{Bf)^o + Bifi - 


f ^2,0), 




B2,0 


— xBofi{Bi^o + Bi^i - 


f ^2,1), 




Bis 


— -Bq^o + xBifi{Bofi - 


f Si.o4 


- B2.0), 


B2S 


= Bi^o + xBi^o{Bifi - 


f Si.i 4 


-B2.1), 


B2,2 


^ Bi^i + xB2,oiBi.o - 


f Bi.i 4 


-B2,l). 



Solving the above system of equations we obtain that the generating function i?o.o — -^(i,i,i)(2;) for the 
non-colored Motzkin numbers of rank 3 satisfies 

= 1 - (1 + x)Bofl + SxS^ o + x^{l- 2x)Bl„ + 2x^5^ q " + x)bIo + x^B^o- 

Now, we would like to give a recursion relation for the sequence ml2'^''^\ This can be automatically done 
with maple's program gfun; The procedure "algeqtodiffeq" translates the (algebraic) equation for _Bo.o = 
M(^i^i^i){x) into an equivalent differential equation of order 7 (it is too long to present here) and then the 
procedure "diffeqtorec " translates the differential equation into a 28-term recursion relation. The first few 
values ofm^n'^''''^ are given in Table\^ this sequence seems not to be listed in 

Example 3.6. Theorem \3.2\ with (|lip gives for r — 4 a set of ^ = 10 equations for the Bi j. As above, it 
is possible to solve this system and obtain that the generating function Bq q = x,i)(2;) for the non-colored 
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71 = 1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


(1,1.1) 

rUn ■ 


1 


4 


16 


78 


404 


2208 


12492 


72589 


430569 


2596471 



Tabic 2: The first few values of the non-colored Motzkin numbers of rank three. 



Motzkin numbers of rank 4 satisfies the following equation 

= l + ix-l)Bo,o-2xBla^x{x + 2){x-l)Blo-x'^{x-2){x + 2)Blf^ 

W{x - l)Blo + x^{x - 2){x + l)^Bo,o + ^^(^ + - 1)^^0,0 " x^'^x^ - 3x - 4)B^o 
+x^{x + l){x - if Bio + x\x - 2){x + XfB^^o + x^{x - l)Bl]o - x^°{x - 2){x + 2)BI% 
-x^\x + 2){x - l)Bl% - 2x^^Bl% + x^'^ix - l)Bl% + x^^BI%. 

The first few values of the corresponding Motzkin numbers ra^'^'^^ of rank 4 O'fe given in Table \3i this 
sequence seems not to be listed in JTI. 





n = 1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


(1.1.1) 


1 


5 


25 


155 


1025 


7167 


51945 


387000 


2944860 


22791189 



Table 3: The first few values of the non-colored Motzkin numbers of rank four. 



Remark. Let us consider the non-colored Motzkin numbers of rank r (i.e., all weights are equal to 1). They 
satisfy a r(r)-term recursion relation where T(r) is defined by this property (and is minimal). This yields a 
well-defined sequence {r(r)}rGN starting - according to the above examples - with 3, 7, 28, . . .. The original 
conjecture T(r) = 2r + 1 proved to be too naive, but it seems that an exact formula for r(r) will be very 
difficult to obtain. Thus, approximations or bounds would also be interesting. Considering the equations for 
the generating functions -Bo,o = .^(i,i,i)(2;) of non-colored Motzkin numbers given in the above examples for 
rank r < 4, one observes that they all have the form /(x,i3o.o) = where f{x,y) is a polynomial of degree 
2^ in y with coefficients in Z[a;] (and with constant coefficient 1); more precisely, the coefficient ai{x) of 
is a polynomial of degree at most i in x over the integers Z. Thus, we can write 

2'^ i 

0=l + ^^a,,,a;^"(Bo.or (12) 

1=1 j=0 

where a^j- € Z. It would be interesting to find out whether this representation holds for all ranks (which is 
what we expect) or whether it is confined to r < 4. 

We would like to close this paper by stressing that it is still an open problem to derive a recursion relation 
for the Motzkin numbers of higher rank. As the explicit examples r — 2,3,4 and the complicated set of 
equations for the generating function given in Theorem 13 . 21 show this will be a rather daunting task (see also 
the last remark). 
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